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Abstract 

A detailed simulation of air showers produced by dust grains has been 
performed by means of the aires Monte Carlo code with the aim of comparing 
with experimental data. Our analysis indicates that extensive dust grain 
air showers must yet be regarded as highly speculative but they cannot be 
completely ruled out. 
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It has long been known that small solid particles (or dust grains) may be accelerated 
effectively at strong collisionless shock waves - popularly in supernovae - In addition, 
it was suggested years later that, even upon destruction, a substantial fraction of the super- 
thermal debris (those created in the pre-shock region) might re-enter the shock waves to 
undergo additional acceleration [g]. These ideas have been recently strengthened by astro- 
nomical observations [H and by the analysis of primitive meteorites @]. 

Despite some still open questions concerning the timescale for grain destruction in the 
source environment and, later on, during the trip to Earth, from time to time dust grains 
have been considered as possible progenitors of giant air showers. The earliest significant 
contribution we are aware of is due to Herlofson . He argued that relativistic dust grains 
could indeed produce extensive air showers, provided that in the first interaction each con- 
stituent nucleon contributes with an energy above 10^^ eV. The field then lay fallow for 
eighteen years until Hayakawa tried to explain the exceptional event reported by Suga et 
al. 1^. Prompted by this proposal several plausible suggestions came out in the 1970s [|]. 

Unfortunately, the values of the depth of shower maximum registered by Haverah Park 
and Volcano Ranch experiments do not provide an exact picture of showers initiated by 
dust grains 0. Nonetheless, in view of the low statistics at the end of the spectrum and 
the wide variety of uncertainties in these experiments, one may be excused for reserving the 
judgment. In order to increase the statistics significantly, the Southern Auger Observatory 
is currently under costruction: A surface array (that will record the lateral and temporal 
distribution of shower particles) + an optical air fluorescence detector (which will observe 
the air shower development in the atmosphere) [0. A major advantage of the optical device 
is precisely its capability of measuring the depth for the maximum shower development. 

Whatever the source(s) of the highest energy cosmic rays, the end of the spectrum 
remains unexplained by a unique consistent model |Tl|] . This may be due to the present lack 
of precision of our knowledge of these rare events, or perhaps, may imply that the origin 
and nature of ultrahigh energy cosmic rays have more than one explanation. If one assumes 
the latter hypothesis, relativistic specks of dust are likely to generate some of the events. 
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The above considerations have motivated us to re-examine the effects of giant dust grain 
air showers. Before proceeding to discuss in detail how these extensive showers evolve, 
it is useful to review some key characteristics of the mechanisms that could lead to the 
destruction of dust grains. Notice that the opacity of the interstellar medium will impose 
lower and upper limits on the value of the Lorentz factor 7 for the primaries impacting on 
the Earth atmosphere. 

An early investigation by Berezinsky and Prilutsky indicates that the grains turn out to 
be unstable with respect to the development of a fracture |]12|. On the one hand, subrela- 



tivistic dust grains disintegrate in collissions with heavy nuclei of the interstellar gas. On 
the other hand, electrical stress induced by the photoelectric effect in the light field of the 
Sun results in the mechanical disruption and subdivision of relativistic grains. Doubts have 
also been expressed about the prospects of surviving against heating arising from photoion- 
ization within the solar radiation field [0. The evaporation of the surface atoms |jT^ and 



even the capture of electrons from the interstellar medium have been suggested as possible 
ways to reduce the accumulation of charges. 

All in all, the path length up to the first break-up in favor of these figures, (log 7 ~ 2 and 
initial radii between 300 - 600 A) turns out to be of a few parsecs |T^, i.e., much less than 



the characteristic size of the Milky Way. This entails that only RX J0852. 0-4622, the closet 
young supernova remnant to Earth (distance ~ 200 pc) |jl6|], could be considered scarcely 
far away. 

Let us now turn to the discussion of dust grain air showers (DGASs). Relativistic dust 
grains encountering the atmosphere will produce a composite nuclear cascade. Strictly 
speaking, each grain evaporates at an altitude of about 100 km and forms a shower of 
nuclei which in turn produces many small showers spreading over a radius of several tens 
of meters, whose superposition is observed as an extensive air shower. For 7^1, the 
internal forces between the atoms will be negligible. What is more, the nucleons in each 
incident nucleus will interact almost independently. Consequently, a shower produced by a 
dust grain containing n nucleons may be modelled by the collection of n nucleon showers. 



each with 1/n of the grain energy. Thus, recalhng that muon production in proton showers 
increases with energy as E^'^^ |]T^, the total number of muons produced by the superposition 



of n individual nucleons showers is, A^^*"^ oc n{E/n)^'^^ , or, comparing to proton showers, 
^^'^ = n^'^^N'P. Of course, these estimates are approximations, but without cumbersome 
numerical calculation one could naively select the event recorded at the Yakutsk array (May 
7, 1989) as the best giant DGAS candidate in the whole cosmic ray data sample [|l^ , [l9|j . 

In order to go a step further and test these qualitative considerations, we have performed 
several atmospheric cascade development simulations by means of the aires Monte Carlo 
code [^]. Several sets of showers were generated, each one for different specks of metallic 



nature, i.e. with different Loretz factors. The sample was distributed in the energy range of 
10^^ eV up to 10^° eV and was equally spread in the interval of 0° to 60° zenith angle at the 
top of the atmosphere. All shower particles with energies above the following thresholds were 
tracked: 750 keV for gammas, 900 keV for electrons and positrons, 10 MeV for muons, 60 
MeV for mesons and 120 MeV for nucleons and nuclei. The particles were injected at the top 
of the atmosphere (100 km.a.s.l), and the surface detector array was put beneath different 
atmospheric densities selected from the altitude of cosmic ray observatories (see Table I). 
SIBYLL routines ETI] were used to generate hadronic interactions above 200 GeV. Notice that 



while around 14 TeV cm. energies the kinematics and particle production of minijets might 
need further attention, for DGASs the energy of the first interaction is reduced to levels 
where the algorithms of sibyll accurately match experimental data p2[. Results of these 
simulations were processed with the help of aires analysis program. Secondary particles of 
different types and all charged particles in individual showers were sorted according to their 
distance R from the shower axis. 

In an attempt to examine our qualitative considerations we first restrict the attention 
to the highest event reported by the group at Yakutsk. Next, details of the most relevant 
observables of the showers are given from the analysis of both the particles at ground and 
those generated during the evolution of the cascade. 

The Yakutsk experiment determines the shower energy by interpolating and/or extrap- 



dating the measurements of peoo (charged particle density at 600 m from the shower core), 
a single quantity which is known from the shower simulations to correlate well with the 
total energy for all primary particle types. In the case of the event detected on May 1989, 
the trigger of 50 ground-based scintillation detectors at 200 - 2000 m from the shower core 
allowed to estimate a reliable value of peoo = 54 m~^ and a declination axis given by cos 9 = 
0.52 [l^. Examining the lateral distribution of our simulation sample, we note that showers 
initiated by relativistic particles (log 7 = 4- 3.8) of energies in the range 36 - 38 EeV and 
with orientation 6 = 59° are compatible with such a value of peoo HH- It is important to 



stress that the lower and upper energy bounds are compatible with the event of interest 
within 2 a; however, for these high Lorentz factors the interstellar medium is extremely 
opaque to dust grains (again the reader is referred to Ref. |1^). In Fig. 1 we show the 



lateral distributions of muons and charged particles. It is easily seen that the predicted 
fluxes at ground level are partially consistent with those detected at the giant array (see 
Fig. 2). 

In what follows we shall discuss the main properties of DGASs. The atmospheric depth 
-^max at which the shower reaches its maximum number of secondary particles is the standard 
observable to describe the speed of the shower development. In view of the superposition 
principle this quantity is practically independent of n, depending only upon 7. For this 
reason the altitude of the X^ax is generally used to find the most likely mass for each primary. 
The behavior of the Xmax with increasing 7 in vertical showers has been already discussed 
in detail elsewhere 0. Here we shall extend Linsley's analyses studying the dependence 
of the shower evolution with the incident angle. In the last panel of Fig. 3 we show the 
numerical results from 10^^ eV cascades initiated by nucleons with log 7 = 3 and different 
primary zenith angles. It can be observed that, at the same total energy, an air shower 
induced by particles with an oblique incidence develops faster than a vertical shower. Since 
muons are typically leading particles in the cascade, the position of Xmax is also related 
to the relative portion of muons in the shower. To illustrate this last point, the resulting 
fluxes at ground level from the same events are shown in first two panels of Fig. 3. As 



it is expected, the radial distribution of the shower particles of inclined primaries (mainly 
dominated by muons) is flatter than the distribution of a vertical shower. The opposite 
behavior points onto the supremacy of electrons and positrons near the core. The density 
of these charged particles, however, falls off rapidly with increasing core distance, dimming 
the electromagnetic cascade. 

To analyse the sensitivity of ground arrays to the shower parameters, the radial variation 
of different groups of secondary particles (we have considered separately 7, e^e~, and /U^/i") 
was studied at two observation levels. In Fig. 4 we show the last steps of the evolution of 
the lateral distribution along the longitudinal shower path. It can be seen that there is 
practically no change in the radial variation. Thus, we conclude that the flux of particles 
does not have intrinsic sensitivity to the observation altitude until 1400 m. 

Coming back to the general features of the longitudinal development, the last exercise is 
to analyse the dependence of the atmospheric shower profile with 7. To this end we carry 
out numerical simulations of vertical showers induced by primaries with the same mass and 
different Lorentz factors, namely log 7 = 2.5 to 3.2. The results are shown in the last 
panel of Fig. 4. As expected, for the same primary mass the number of secondary particles 
increases with rising 7, but it is interesting to note that both cosmic ray cascades present 
similar shapes and peak around the same atmospheric depth, i.e., Xmax ~ 350 ± 47 g/cw? 
(consistent with the analysis of 0). 

Putting all this together one can draw the following tentative conclusions: 
(i) The Fly's Eye collaboration has presented evidence indicating that typical extensive air 
showers above ~ 1 EeV develop at a rate which is consistent with a steep power law spectrum 
of heavy nuclei that is overtaken at higher energies by a flatter spectrum of protons . The 
group of the Akeno Giant Air Shower Array has reported 7 events above 10^° eV until August 



1998 p5|. In general, the muon component agrees with the expectation extrapolated from 



lower energies |2^. However, the highest energy event recorded at Yakutsk (E' = 1.1 ± 0.4 
eV) seems to be a rare exception. It can be readily noticed from Fig. 1 that the almost 
completely muonic nature of particles detected may be explained by a 36 - 38 EeV DGAS. 



Besides, if this speculation is true the lack of obvious candidate sources at close distances 
leaves the origin of this event still as a mystery. 

(ii) The dependence of the DGAS longitudinal profile on 7 in the studied range is rather 
weak. The shower disc becomes slightly flatter and thicker with decreasing primary energy 
and/or increasing zenith angle. This is mainly due to the rising content of muons among all 
charged particles in the inclined shower. 

Dust is a very widespread component of diffuse matter in the Galaxy with apparently 
active participation in charge particle acceleration. If such acceleration is possible, dust 
grains may play an important role in determining the composition of galactic cosmic rays. 
The question of whether the opacity of the interstellar medium could prevent relativistic 
dust grains from reaching the Earth is as yet undecided. Observation of giant DGASs would 
give an experimental and definitive answer to this question. We recommend that Auger 
search data be analysed for evidence of DGASs. 

I am grateful to D. Gomez Dumm for helpful comments on the original manuscript and 
to F. Zyserman for computational aid. I am also indebted to Yakutsk Collaboration, and 
particularly to M. I. Pravdin for permission to reproduce Fig. 2. This work has been 
partially supported by CONICET. The author is on leave of absence from UNLP. 
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TABLES 



TABLE I. Sites of giant air shower arrays 



Experiment 



longitude 



latitude 



altitude 



Volcano Ranch 

SUGAR 

Haverah Park 

Yakutsk 

Fly's Eye 

AGASA 

Auger 



35°09'N 

30°32'S 

53°58'N 

61°42'N 

40° 'N 

38°47'N 

35°12'S 



106°47'W 

149°43'E 

1°38'W 

129°24'E 

113°'W 

138°30'E 

69°12'W 



834 g/cm2 
1020 g/cm^ 
1020 g/cm^ 
1020 g/cm^ 
860 g/cm^ 
920 g/cm^ 
875 g/cm^ 
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FIG. 1. From left to right, (a) Lateral distributions of charged particles and muons from AIRES 
simulations of a 36 EeV, log 7 = 4 DGAS. The error bars indicate the RMS fluctuation of the 
means, (b) Id. with E = 38 EeV and log 7 = 3.8. 
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FIG. 2. Density of muons recorded at Yakutsk by ground and underground detectors on May 7, 
1989 at 13 h 23 min Greenwhich time. There is a remarkable contradiction with an extrapolation 
from results of the low energy region (dashed line). This figure was originally published in Ref. 
[18]. 
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FIG. 3. Prom left to right: (a) Lateral distribution of charged particles of DGASs of 10 EeV. 
The primaries with log 7 = 3 were injected vertically and with 60 degree zenith angle. The 
figure also shows the fluctuations measured in terms of RMS. (b) Idem for the case of muons. (c) 
Longitudinal evolution of the shower profiles. The error bars indicate the standard fluctuations. 
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FIG. 4. Lateral distributions of electrons, muons and gammas (at ground level) from AIRES 
simulations of 10^^ eV vertical showers (7 = 300). The figure also shows the longitudinal cascade 
development profile of 10^^ - 10^^ eV showers as would be seen by the Auger fiuorescence detector. 
The logarithm of the Lorentz factors are 2.5 and 3.2 respectively. The error bars indicate the 
standard fluctuations (the RMS fluctuations of the means are always smaller than the symbols) . 
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